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ABSTRACT 


Results  of  an  experimental  study  on  the  dynamic  properties 
of  6061-T6  aluminum  are  presented.  Areas  studied  included 
stress-strain-strain  rate  and  reverse  loading  behavior,  elas¬ 
tic  constants,  equation  of  state,  compressive  and  release 
wave  characteristics,  and  spall  fracture.  The  material  showed 
approximately  elastic-perfectly  plastic  behavior  under  uni¬ 
axial  stress  compression  (0.005  to  1000/sec  strain  rates) . 
Strain- rate  sensitivity  was  very  low  at  20 °C,  but  increased 
at  temperatures  of  150,  204,  and  288°C.  Longitudinal  and 
shear  wave  velocities  at  20°C  were  6.368  and  3.197  mm/ysec, 
and  temperature  and  pressure  dependence  was  also  measured  and 
various  elastic  constants  were  calculated.  The  shock  wave 

equation  of  state  (EOS)  up  to  200  kbar  was  determined  as 

2 

au  =  1.0  +  140,4  u  +  37.7  u  .  The  EOS  was  also  calculated 
H  P  P 

from  the  elastic  constants  and  comparisons  made  with  the  shock 
wave  EOS.  Compressive  wave  tests  showed  a  well-defined  elastic 
precursor  with  a  steady-state  value  of  ^  5  kbar.  The  yield 
level  from  shock  wave  profile  measurements  was  approximately 
the  same  as  that  determined  from  uniaxial  stress  and  EOS  tests. 
The  impact  velocity  required  for  spall  fracture  was  found  to 
increase  with  decreasing  pulse  width,  and  to  decrease  with  in¬ 
creasing  temperature. 
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FOREWORD 

The  prediction  of  reentry  vehicle  response  to  impulsive  loading 
resulting  from  energy  deposition  has  been  studied  extensively 
during  the  last  decade.  Analytical  models  and  computer  routines 
have  been  developed  to  assess  the  vulnerability  of  such  vehicles 
to  the  initial  loading  phase,  in  which  stress  waves  are  gener¬ 
ated  and  propagated  through  the  structure,  and  to  subsequent 
phases,  where  elastic  vibration,  plastic  deformation  and/or  frac¬ 
ture  may  occur. 

A  Defense  Atomic  Support  Agency  (DASA)  program  was  initiated 
with  the  objective  of  improving  computer  codes  for  the  predic¬ 
tion  of  damage  induced  by  X-rays  (PREDIX) .  The  PREDIX  metals 
program  has  combined  the  efforts  of  several  contractors  under 
the  direction  of  DASA.  The  primary  contributions  of  each  con¬ 
tractor  can  be  summarized  as  follows: 


General  Motors  Corporation — Measurement  of  material 
properties,  including  stress-strain-strain  rate  be¬ 
havior,  elastic  constants,  equations  of  state,  shock 
wave  profiles  and  spall  fracture. 

Effects  Technology  Inc. --Spall  fracture  tests  and 
analysis  (exploding  foil)  and  underground  testing. 

KMS  Technology  Center--Cons titutive  modeling. 

Physics  International  Company — Energy  deposition  and 
spall  fracture  studies  (electron  beam) . 

Systems,  Science  and  Sof tware--Energy  deposition  and 
front  surface  phenomena,  constitutive  modeling  and 
code  development,  and  underground  testing. 
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INTRODUCTION 


This  report  is  one  of  a  series  of  six  giving  results  of  a  study 
of  dynamic  properties  of  four  metals.  The  principle  objectives 
of  this  study  were: 

•  Provide  measurements  of  dynamic  properties  of 
materials  to  be  used  as  inputs  to  model  develop¬ 
ment  for  shock  wave  propagation  and  dynamic 
fracture. 

•  Provide  data  to  establish  accuracy  of  computer 
code  predictions. 

•  Provide  physical  interpretation  of  experimentally 
observed  material  response  to  dynamic  test  condi¬ 
tions  . 

The  primary  areas  studied  included:  (1)  Stress-strain-strain 
rate  behavior  in  uniaxial  stress,  including  Bauschinger  effect; 
(2)  Elastic  constants,  including  temperature  and  pressure 
dependence;  (3)  Hugoniot  equation  of  state;  (4)  Compressive 
and  release  wave  profiles,  including  elastic  precursor  decay 
and  wave  attenuation;  and  (5)  Spall  fracture  studies,  including 
pulse  width  dependence.  The  six  reports  in  this  series  are: 


Vol. 

I: 

Summary  of  Results,  DASA  2501-1  (AD712847) 

vol. 

II: 

Experimental  Methods  and  Techniques, 

DASA  2501-2 

Vol. 

III: 

6061-T6  Aluminum,  DASA  2501-3 

Vol. 

IV: 

Alpha  Titanium,  DASA  2501-4 

Vol. 

V: 

OFHC  Copper,  DASA  2501-5 

Vol. 

VI: 

Tantalum,  DASA  2501-6 
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The  experimental  techniques  utilized  in  this  study  are  reviewed 
in  Vol.  II  (DASA  2501-2) .  Briefly,  the  principal  techniques 
applied  were: 

Stress-strain  Studies — A  laboratory-type  universal 
testing  machine  was  used  for  rates  <  0.1/sec,  a 

_3 

medium  strain  rate  machine  for  rates  of  10  /sec 
2 

to  10  /sec,  and  a  Hopkinson  bar  device  for  rates 
2  3 

of  about  10  to  5  X  10  /sec. 

Elastic  Constants — The  pulse  superposition  method  was 
used  for  measuring  longitudinal  and  shear  wave  velo¬ 
cities  as  functions  of  temperature  and  hydrostatic 
pressure. 

Equation  of  State — Gun-Launched,  flat-plate  impact 
techniques  were  used  to  generate  uniaxial  strain 
conditions.  Velocities  up  to  0.6  mm/ysec  were 
achieved  with  63.5  mm  and  102  mm  single-stage  com¬ 
pressed-gas  guns,  and  velocities  of  0.3  to  2  mm/ysec 
were  achieved  with  a  64  mm  single-stage,  powder  gun. 
Hugoniot  data  were  obtained  with  x-cut  quartz  gage 
and  streak  camera  techniques. 

Wave  Propagation  and  Spall  Fracture — Flat-plate  impact 
techniques  were  also  used.  Wave  profiles  were  measured 
with  x-cut  quartz  gages  and  a  laser  velocity  interfero¬ 
meter. 

A  literature  search  on  dynamic  properties  of  materials  was 
carried  out  under  the  PREDIX  program  and  a  listing  of  recent 
publications  (including  abstracts)  on  dynamic  properties  of 
materials  is  given  in  Reference  1. 
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SECTION  I 

MATERIAL  PROPERTIES 


Aluminum  has  a  face- centered- cubic  structure  and  the  6061  type 
on  which  the  tests  discussed  in  this  report  were  carried  out 
is  a  heat-treatable  wrought  alloy.  This  material  has  nominally 
3%  alloying  constituents  (primarily  magnesium,  silicon  and  iron) 
and  results  of  chemical  analysis  are  given  in  Table  I.  The 
alloy  is  solution-heat- treated  and  precipitation-hardened  to 
the  ' -T6'  temper,  as  outlined  in  Table  I.  The  material  was 
purchased  as  sheet  and  plate  stock  in  thicknesses  of  1.6,  3.2, 
6.4,  12.7  and  50.8  mm.  Although  all  material  was  ordered  to 
the  same  specifications,  there  were  slight  differences  in  yield 
and  flow  stress  behavior  for  the  different  thicknesses.  Hard¬ 
ness  of  the  various  thicknesses  was: 


Thickness 

Measured 

Estimated 

(mm) 

rb 

Bhn  (500  kg) 

1.6 

%  50 

83 

3.2 

51 

84 

6.4 

57 

91 

12.7 

59 

94 

50.8 

59 

94 

All  stock  thicknesses  howed  the  same  degree  of  grain  structur¬ 
ing,  although  there  was  some  variation  in  grain  size  with  the 
thinner  stock  having  smaller  grains.  Photomicrographs  showing 
the  grain  structure  in  the  3.2  mm  and  12.7  mm  stock  are  shown 
in  Figure  1. 
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TABLE  I 

6061  ALUMINUM  CHARACTERISTICS 


COMPOSITION  (WEIGHT  %) : 


Element 

Maximum 

Measured 

(Ref.  2) 

(12.7  mm  Plate) 

Silicon 

0.80 

0.63 

Iron 

0.70 

0.41 

Copper 

0.40 

0.20 

Manganese 

0.15 

0.06 

Magnesium 

1.20 

0.84 

Zinc 

0.25 

0.09 

Chromium 

0.35 

0.20 

Titanium 

0.15 

0.06 

PROCESSING  (REF.  2): 

A.  SOLUTION  HEAT  TREATMENT,  -T4  TEMPER 

1.  970°F  for  0.2-1  hours  in  molten  salt. 

2.  Quench  in  cold  water  with  minimum 
delay. 


B.  PRECIPITATION  HEAT  TREATMENT,  -T6  TEMPER 

1.  350°F  for  6-10  hours. 

2.  Cooling  rate  unimportant  (not  too  slow). 


C.  STRESS -RELIEVED  BY  STRETCHING,  -T651  TEMPER 
(OPTIONAL) 

1.  1  1/2-3%  permanent  set  after  B. 
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3 

The  average  measured  density  was  2.703  g/cc  (0.0976  lb/in  ). 

For  use  in  equation  of  state  and  energy  deposition  calculations, 
several  physical  constants  were  compiled  from  the  li terature 
and  are  listed  below: 

Volume  coefficient  of  thermal  expansion,  3: 

67.0  x  10~6  +  0.06  x  10~6  T  /°C 

Specific  heat,  c  : 

P 

0.220  +  0.00015  x  T  cal/g-°C 
Solidus  temperature:  616°C 
Liquidus  temperature:  651°C 
Boiling  temperature:  ^  2450°C 
Latent  heat  of  fusion:  'v  9  3  cal/g 
Latent  he  t  of  vaporization:  2510  cal/g 

(T  in  °C) 

No  change  in  the  microstructure  of  6061-T6  aluminum  was  ob¬ 
served  for  shock  stresses  up  to  10  kbar,  i.e.,  there  was 
no  visual  evidence  of  twinning  or  grain  deformation.  Photo¬ 
micrographs  of  as-received  and  shocked  material  are  compared 
in  Figure  2. 


6 


( 

F 

S 


t 

! 


MANUfACTURlNG  DEVELOPMENT 


•  GENERAL  MOTORS  CORPORATION 


MSL-70-23,  Vol.J'I 


TJ 

c 

(0 


n 

> 


•ri 

<u 

o 

<0 

K 

I 

tn 

< 


E  -i 

3 

C  H 
•H  M 

E  0) 

3  4J 

■H  (0 

<  £ 


II 


VO  -O 
Eh  <L) 

i  .* 
-t  u 
vo  o 

O  jC 

vo  tn 


rv| 

ffl 

lH 

3 

0> 

*H 

Cn 


7 


L 


A 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 


MSL-70-23,  Vol.  HI 


SECTION  II 

STRESS-STRAIN  STUDIES 


Shock  wave  and  structural  response  code  input  requirements  in¬ 
clude  constitutive  equations  relating  stress  to  strain,  strain 
rate  and  temperature.  To  obtain  these,  measurements  are  made 
of  uniaxial  compressive  stress  vs.  strain  at  various  strain 
rates  and  temperatures.  These  tests  give  yield  and  flow 
stress  behavior,  strain  rate  sensitivity  and  work  hardening 
characteristics.  Although  it  is  recognized  that  heating  rate 
and  time-at-temperature ,  as  well  as  final  temperature,  are 
important  in  assessing  material  properties  under  short-time 
(<  10  ®  sec)  energy  deposition  conditions,  most  stress-strain 
data  reported  here  are  for  relatively  long  heating  times  only. 
Bauchinger  effect  tests  were  conducted  to  provide  information 
on  unloading  and  subsequent  yield. 


UNIAXIAL  STRESS  TESTS 


Results  of  compressive  stress-strain  tests  on  6061-T6  are 
shown  in  Figures  3  and  4  for  strain  rates  of  0.005/sec  to 
10  /sec  and  test  temperatures  of  20°  to  288cC,  with  each 
curve  the  average  of  several  tests.  At  20°C,  strain  rate 
sensitivity  is  negligible  (within  normal  material  scatter 
and  test  errors) ,  but  increasing  the  temperature  above  150 °C 
results  in  an  increase  in  rate  sensitivity,  yield  stress  and 
work  hardening  decrease  with  increasing  temperature. 
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Figure  4  Compressive  Yield  Stress  vs.  Log 
Strain  Rate,  6061 
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The  temperature  data  in  Figures  3  and  4  were  obtained  for  a 

heating  rate  of  ^  0.5°C/sec  and  it  has  been  found  that  yield 

and  flow  stress  at  elevated  temperatures  are  a  function  of 

( 8 ) 

heating  rate  and  time-at-temperature .  '  For  example,  the 

tensile  yield  stress  in  6061-T651  is  shown  in  Figure  5  as  a 
function  of  heating  rate  for  two  final  temperatures.  The  data 
indicate  that  degradation  of  yield  at  elevated  temperature  has 
three  distinct  regimes.  At  high  heating  rates  (>  l°C/sec  for 
260°C)  the  yield  shows  little  dependence  on  heating  rate,  sug¬ 
gesting  an  almost  instantaneous  decrease  in  yield  as  tempera- 

-4 

ture  increases.  At  intermediate  rates  (10  to  l°C/sec)  a 

time-dependent  softening  occurs  and  the  yield  decreases  with 

decreasing  heating  rate  due  to  overageing  or  precipitate  coar- 

-4 

sening.  At  low  heating  rates  (<  10  °C/sec)  the  yield  reaches 
a  constant  value.  Although  the  data  in  Figure  5  are  presented 
in  terms  of  heating  rate,  the  conclusions  would  be  qualitatively 
the  same  if  discussed  as  a  function  of  time-at-temperature. 

All  tests  were  performed  at  constant  strain  rate  and  the  ob¬ 
served  response  may  be  rate  dependent. 

There  was  a  slight  variation  in  yield  and  flow  stress  for  test 
specimens  taken  from  stock  of  different  thicknesses,  as  shown 
in  Figure  6.  This  variation  is  comparable  to  that  due  to 
scatter  from  one  plate  alone  and  to  rate  sensitivity  at  20°C. 

REVERSE  LOADING  TESTS 

★ 

The  reverse  loading  behavior  or  Bauschinger  effect  was  studied 
by  performing  uniaxial  tensile  stress  tests  after  the  material 
had  been  prestrained  in  the  opposite  (compressive)  sense. 

TE - 

The  Bauschinger  effect  can  be  defined  as  the  reduction  in  stress 
at  a  given  strain  in  the  reverse  direction  compared  to  that  in 
the  original  loading  direction. 
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Figure  6  Compressive  Stress  vs.  Strain, 
Material  Effect,  6061-T6 
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Specimens  were  loaded  in  compression  to  strains  of  0.5  to  3.0%, 
unloaded  to  zero  stress,  and  then  immediately  loaded  in  ten¬ 
sion.  The  Bauschinger  effect  can  be  evaluated  in  terms  of  the 
ratio  of  yield  stress  in  tension  (reload)  to  the  yield  stress 
in  compression  or  to  the  maximum  preload  stress.  Care  must  be 
taken  in  interpreting  data  from  tests  where  small  (<  0.2%) 
values  of  strain  are  measured.  Choosing  a  definition  such  as 
the  first  deviation  from  a  linear  elastic  behavior,  or  stress 
at  a  given  offset  strain,  will  result  in  various  interpretations 
of  the  Bauschinger  effect.  Reverse  loading  test  results  are 
shown  in  Figure  7  and  the  data  is  cross-plotted  in  Figure  8 
where  two  definitions  of  the  Bauschinger  effect  are  plotted 
against  plastic  prestrain. 

(9) 

Following  Milligan,  Koo  and  Davidson  ,  a  Bauschinger  effect 
factor  (BEF)  can  be  defined  in  two  ways: 

(1)  The  ratio  of  yield  stress  upon  reverse  loading 
to  the  initial  yield  stress  (see  Figure  8  for 
symbol  definition) 

°C 

BEF  =  — 

°A 

(2)  The  ratio  of  yield  stress  during  reloading  to 
the  maximum  stress  upon  loading 

o 

BEF  ■  — 

°B 


This  latter  definition  takes  into  account  strain  hardening. 
All  stresses  were  defined  at  a  0.1%  offset  strain.  From 
Figure  8  it  appears  that  the  Bauschinger  effect  factor  de¬ 
creases  (Bauschinger  effect  increases)  up  to  about 
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Figure  8  Bauschinger  Effect  Data,  6061-T6 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 


MSL-70-23,  Vol.  Ill 

1.0  -  1.5%  plastic  strain.  Then,  with  increasing  prestrain, 

the  BEF  from  both  definitions  appears  to  increase  in  magnitude. 

These  results  are  quantitatively  similar  to  the  results  on 

6061-T6  plate-stock  aluminum  and  modified  high-strength 
(9 ) 

steel  reported  by  Milligan,  et  al.,  although  the  preload 
was  in  tension  rather  than  compression.  The  probable  cause 
for  the  increase  in  BEF  after  reaching  a  minimum  at  low  strains 
has  been  attributed  to  strain-hardening  ef fects .  Fig¬ 

ure  9  shows  results  from  a  set  of  cyclic  compression- tens  ion 
tests  on  6061-T6  aluminum.  The  primary  conclusion  to  be  drawn 
from  this  data  is  the  type  of  hardening  mechanism  for  this  mate¬ 
rial.  At  a  strain  of  approximately  0.5%,  kinematic  hardening 
predominates.  ^ 

BIAXIAL  STRESS  TESTS 

Tubular  specimens  machined  from  6061-T6  bar  stock  were  tested 
under  biaxial  stress  conditions  to  determine  the  yield  locus 
(or  surface) .  A  0.25%  strain  offset  technique  in  the  octahedral 
shear  stress-strain  plane  was  used  to  determine  yield.  Results 
are  given  in  Figure  10  and  show  that  both  the  von  Mises  and  the 
Tresca  yield  theories  provide  reasonable  fits  to  the  data. 


(Stress-strain  tests  of  6061-T6  have  been  conducted  under  a 
number  of  other  programs  at  General  Motors  Corporation.  The 
reader  is  referred  to  previously  published  reports  and  papers 
for  additional  details  in  such  areas  as  uniaxial  stress  test¬ 


ing 


(12-14)  biaxial  stress  testing (8 , 15 , 16) t  heating  rate 


effects ,  and  preshocking  and  prestraining  ef f ccts (19) . ) 
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Figure  9  Cyclic  Loading  Behavior,  6061-T6 
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SECTION  III 

ELASTIC  BEHAVIOR 


Measurement  of  elastic  constants  provides  additional  inputs  to 
the  study  of  material  response.  The  three  basic  measurements 
are  density,  longitudinal  (dilatational)  wave  velocity  and 
shear  (transverse)  wave  velocity.  Values  obtained  at  standard 
temperature  and  pressure  (20°C  and  atmospheric  pressure)  estab¬ 
lish  the  parameters  used  in  first-order  (temperature  and  pres¬ 
sure  independent)  calculations  of  material  response.  Measure¬ 
ments  of  temperature  and  pressure  dependence  permit  more  exact 
calculations  to  be  made  and,  in  particular,  lead  tc  prediction 
of  pressure-compression  isotherms,  adiabats  and  hydrostats.  If 
these  measurements  are  made  with  sufficient  accuracy  (<  0.1%), 
the  calculated  pressure-compression  behavior  can  frequently 
be  extrapolated  with  reasonable  confidence  to  higher  pressures 
than  covered  in  the  actual  measurements. 


WAVE  VELOCITY  MEASUREMENTS 

Initial  measurements  of  longitudinal  and  shear  wave  velocities 
were  made  at  20°C  and  atmospheric  pressure  (i.e.,  P  =  0) . 
Values  obtained  were: 

CT  =  6.368  ±  0.006  mm/usec 

•U 

Cc  =  3.197  ±  0.003  mm/ysec 

O 
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Although  the  accuracy  of  these  measurements  is  estimated  to 
be  about  ±  0.1%,  it  should  be  noted  that  Cr  and  Cc  may  vary 
significantly  depending  on  the  structural  form  of  the  mate¬ 
rial  (i.e.,  plate,  bar,  etc.),  thermal  and  mechanical  history, 
and  the  measurement  direction  relative  to  any  structural  ani¬ 
sotropy.  A  review  of  the  literature  on  elastic  wave  velo¬ 
cities  in  6061-T6  gave  values  of  CT  ranging  from  6.34  to 
6.48  mm/ysec  and  Cg  ranging  from  3.04  to  3.26  mm/ysec. 

These  variations  may  be  partly  due  to  differences  in  experi¬ 
mental  technique  but  seem  to  be  primarily  the  result  of  differ¬ 
ences  in  composition,  structure,  texture,  and  orientation. 

For  example,  with  50.8  mm  plate  stock,  Cr  normal  to  the  plate 
was  about  0.5%  less  than  CL  in  the  plane  of  th  plate. 

Measurements  were  also  made  of  the  temperature  and  hydrostatic 
pressure  dependence  of  the-  wave  velocities.  Correction  of  the 
data  for  thermal  expansion  and  hydrostatic  compression  gave  the 
following: 

Temperature  Dependence,  0-100°C,  P  =  0 
CT  =  6.387  -  (0.00093  T)  mm/ysec 

2j 

Cg  =  3.214  -  (0.00084  T)  mm/ysec 
(T  in  °C) 


Pressure  Dependence,  0-9  kbar,  T  =  20°C 


CT  =  6.368  +  (0.0180  P) 

h 

Cg  -  3.19?  +  (0.0108  P) 
(P  in  kbar) 


2 

(0.0001  P  )  mm/ysec 
2 

(0.00005  P  )  mm/ysec 
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(21) 

For  comparison,  Asay,  et  al.  ,  reported  the  temperature 
dependence  of  CT  and  Cc  in  6061-T6  as  CT  =  6.428  - 
0.00091  x  T  and  Cg  =  3.199  -  0.00086  x  T.  They  also  deter¬ 
mined  the  third-order  elastic  constants  at  25°C  using  a  uni¬ 
axial  stress  technique  and  obtained  =  -3010  kbar, 
v2  =  -1310  kbars  and  vg  =  1410  kbar,  where  v^,  v^,  and 


are  the  third-order  Lame  constants  for  an  isotropic  elastic 

,  .  (24) 

body. 


ELASTIC  CONSTANTS 


* 

Assuming  the  material  is  isotropic  (CL  and  Cg  are  independent 
of  direction  of  measurement  in  a  polycrystalline  specimen) , 
wave  velocity  data  are  used  to  calculate  various  elastic  con¬ 
stants  (see,  e.g.,  Refs.  28-31).  At  20°C  and  zero  pressure, 
the  following  adiabatic  constants  were  obtained  (see  Appendix  A) 


Bulk  wave  velocity,  Cfi  =  5.19  mm/ysec 
Sound  wave  velocity,  Cg  =  5.22  mm/ysec 
Rayleigh  wave  velocity,  C_  =  2.99  mm/ysec 
Poisson's  ratio,  v  =  0.331 
Bulk  modulus,  K  =  728  kbar 
Shear  modulus,  G  =  276  kbar 
Elastic  modulus,  E  =  735  kbar 
Lame*' s  parameter,  \  =  544  kbar 


2C , 


*  ‘”'■'44 

A  cubic  material  is  elastically  isotropic  when 


C  -  C 
^11  12 


=  1- 


where  C^,  C^  and  C12  are  elastic  stiffness  constants.  Single¬ 
crystal  aluminum  (99.95%)  is  only  slightly  anisotropic  with 
2C , 


'44 


C11  '  C12 


^  1.2  at  20°C . 


(25-27) 
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The  bulk  and  shear  moduli  are  used  in  the  development  of  the 
ultrasonic  equation  of  state.  They  were  obtained  first  as 
functions  of  T  and  P  from  the  following: 


This  gave: 

Temperature  dependence,  P  =  0 

K  =  732  -  (0.176  T)  kbar 

Gs  =  279  -  (0.165  T)  kbar 

Pressure  dependence,  T  =  20°C 

Kg  =  728  +  (4.75  P)  -  (0.01  P2)  kbar 

Gg  =  276  +  (2.25  P)  -  (0.006  P2)  kbar 

where  the  superscript.  S  indicates  adiabatic. 

Finally,  isothermal  values  of  K  and  G  as  well  as  adiabatic 
and  isothermal  pressure  and  temperature  derivatives  at  20°C 
and  zero  pressure  were  calculated  using  the  equations  given 
in  Appendix  A.  Results  are  given  in  Table  II.  Use  of  these 
constants  in  calculating  the  isotherm,  adiabat  and  hydrostat 
is  discussed  in  the  section  on  Ultrasonic  Equation  of  State. 

GRUNEISEN  PARAMETER 

The  Gruneisen  ratio  y  is  a  parameter  originally  given  by 
Gruneisen  in  a  solid  equation  of  state  relating  pressure  to 

..  '  23 
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TABLE  II 

ELASTIC  CONSTANTS  FOR  6061-T6  ALUMINUM 
(P  =  0,  T  =  20°C) 


PARAMETER 

VALUE 

KS 

728  kbar 

(3£L)  =  Ks ' 

'3P  ' T  OT 

4.75 

-0.175  kbar/°C 

(H~)S  =  Kos 

4.60 

,|k!, 

l3T  'S 

5.71  kbar/°C 

T 

K 

700  kbar 

(H-)T  =  K0T 

4.87 

'  3T  '  P 

-0.281  kbar/°C 

=  kt’ 

' 3P  ;S  OS 

4.64 

T 

,3k\ 

"5t~  s 

5.75  kbar/®C 

GS  =  GT 

276  kbar 

(3P}T  =  G0T 

2.25 

l 

Wp 

-0.164  kbar/°C 

Ws  “  Gos 

2.12 

(|£j 

2.64  kbar/°C 
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volume  and  energy,  y  =  V  (3P/3E)V  (see,  e.g.,  Refs.  32-34). 

This  parameter  can  be  defined  thermodynamically  as: 


For  an  isotropic  elastic  solid,  this  can  be  expressed  as 

where  CT  and  C0  are  measured  under  adiabatic 

J-l  b 

conditions . 

At  2°°C  and  zero  pressure  this  gave 
yQ  =  1.97 

The  zero-pressure  Gruneisen  parameter  can  also  be  estimated 

(35) 

by  several  other  methods,  including  those  of  Slater  ,  Dug- 

06)  (37  38) 

dale  and  MacDonald  ,  Anderson  and  Dienes  '  '  and  Schuele 

(37  39) 

and  Smith  '  .  These  methods  are  discussed  in  Appendix  B. 

An  estimate  of  the  temperature  dependence  of  the  Gruneisen 
parameter  can  be  obtained  by  differentiating  Equation  3  with 
respect  to  temperature.  At  constant  pressure  (P  =  0,  T  = 

20 °C) : 

»  (st)  p 

(5) 
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The  pressure  dependence  can  be  estimated  by  differentiating 
with  respect  to  pressure  (see  Appendix  C) .  At  constant  temper¬ 
ature  (T  =  20°C,  P  =  0)  : 


•1  t  t 

Another  approach  is  to  assume  y/V  constant/  which  permits 
(3y/3P)T  to  be  estimated  directly  from  the  bulk  modulus  data 
(see  Appendix  C) .  At  constant  temperature  (20°C) : 


(700  +  4.87P)1,205 
=  -0 . 0028/kbar  @  P  =  0 
=  -0.0015/kbar  @  P  =  100  kbar 


DEBYE  TEMPERATURE 


The  Debye  temperature  9  is  a  parameter  related  to  the  maximum 
oscillation  frequency  a  system  assumed  to  be  an  elastic 
continuum  can  support.  It  is  an  important  parameter  in 
thermal  energy  calculations  and  indicates  the  temperature 
above  which  variations  due  to  temperature  for  some  thermodynamxc 
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parameters  such  as  specific  heat,  thermal  expansion  ani  Gru- 
neisen  ratio  become  small.  There  are  a  number  of  methods  for 
calculating  the  Debye  temperature  ^  and  one  approxi¬ 

mation  suitable  for  use  with  elastic  wave  velocity  data  for 
polycrystalline  metals  is: 


where  h  is  Planck's  constant,  k  is  Boltzmann's 
constant,  N  is  number  of  mass  points  and  V  is 
sample  volume. 


At  20°C,  the  elastically-determined  Debye  temperature  for 
6061-T6  is: 


for  N  =  4  (atoms  per  unit  cell),  V  =  (4.04  x 
-7  3 

10  mm)  (lattice  constant),  CT  =  6.3f8  mm/ysec 

Li 

and  Cs  =  3.197  mm/ysec. 

(43) 

Sutton  reported  a  Debye  temperature  for  single-crystal 
aluminum  (99.93%)  of  407°K  at  20°C. 


27 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 


MSL-70-23,  Vol.  Ill 


SECTION  IV 

EQUATION  OF  STATE 


The  development  of  material  response  models  and  computational 
codes  requires  extensive  data  on  material  behavior  under  uni¬ 
axial  strain  conditions,  including  equation  of  state,  wave 
propagation  and  spall  fracture.  Most  data  are  used  directly 
in  developing  models  of  material  behavior,  although  independent 
check  data  (attenuated  wave  shapes  and  spall  wave  profiles)  are 
necessary  to  determine  accuracy  of  the  calculations.  This  sec¬ 
tion  presents  the  equation  of  state  analysis,  and  the  following 
two  sections  cover  wave  propagation  and  spall  fracture. 


SHOCK  WAVE  EQUATION  OF  STATE 

The  hugoniot  equation  of  state  is  the  locus  of  equilibrium 
states  reached  after  shocking  of  a  material,  and  data  are 
usually  obtained  either  as  stress-particle  velocity  points  or 
as  shock  velocity-particle  velocity  poincs  (see,  e.g.,  Refs. 
44-49) .  The  shock  wave  hugoniot  data  for  6061-T6  presented 
in  this  report  were  obtained  with  x-cut  quartz  gages  and  op¬ 
tical  techniques.  Representative  quartz  gage  records  for 
direct  impact  and  transmitted  wave  tests  are  shown  in  Figure  11. 
The  buffered  direct  impact  method  (tungsten  carbide  buffer  plate 
on  front  of  the  quartz)  permitted  stresses  up  to  65  kbar  in  the 
aluminum  while  keeping  the  stress  in  the  quartz  at  an  acceptable 
level.  The  use  of  quartz  gages  for  transmitted  wave  tests  was 
primarily  for  the  study  of  compressive  wave  development  and 
elastic  precursor  decay,  and  results  are  discussed  in  the  Wave 
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DIRECT  IMPACT 
A I  — ►  Q 

V  |  -  0. 139  mm/  /i  s 
^  *  10.8  k b a  r 

DIRECT  IMPACT,  BUFFERED 
AI-^WC/Q 

V  j  -  0.  460  mm/  /is 

O'  -  64.0  kbar 
H 

TRANSMITTED  WAVE 
A I  — A I  /  Q 

V  j  -  0. 173  mm/  /is 
O'  -  13.  4  kbar 


H 


Quartz  Gage  Records 
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Propagation  section.  The  direct  impac*  records,  which  gave  a 
stress-particle  velocity  point  since  the  hugoniot  of  quartz  is 
known,  showed  a  rise-time  of  less  than  10  nsec.  The  slight  round¬ 
ing  at  the  front  of  the  direct  impact  records  in  Figure  11  is  due 
to  impact  tilt,  and  the  ramped  peak  stress  region  results  from 

finite-strain  effects  in  the  quartz  which  are  corrected  in  the 

.  „  .  (50) 

data  analysis . 

Hugoniot  data  at  stresses  of  70  to  160  kbar  were  obtained  by 
using  a  steak  camera  and  a  target  in  the  form  of  a  two-step 
cylinder  ("hat"  configuration) .  A  schematic  of  the  technique 
and  a  sample  record  are  shown  in  Figure  12.  The  two  surfaces 
of  the  target  are  viewed  through  the  camera  slit  and  shock  ar¬ 
rival  is  observed  as  a  change  in  surface  reflectance.  The 
average  shock  velocity  is  obtained  from  the  record,  and  par¬ 
ticle  velocity  is  assumed  to  be  1/2  Vj  (impact  velocity)  since 
impactor  and  target  are  the  same  material. 

The  hugoniot  may  be  expressed  in  several  forms.  A  convenient 
form  for  experimental  work  is  that  established  by  a  least-squares 
stress  and  particle  velocity,  as  established  by  a  least-squares 
fit  tc  the  data  in  the  stress-particle  velocity  (aH  -  up)  plane. 
Transformation  of  the  hugoniot  into  various  planes,  such  as 
shock  velocity-particle  velocity  (Ug  -  u^)  and  stress-volume 
(o  -  v) ,  xs  performed  by  assuming  a  material  model.  This  was 
done  by  assuming  an  ideal  elastic-plastic  wave  structure  with 
equilibrium  initial  and  final  states  and  applying  the  mass  and 
momentum  conservation  equations.  This  leads  to: 
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Therefore,  when  either  au  -  u  or  U„  -  u  relations  are  estab- 

ii  p  bp 

lished  and  a  ,  p  ,  u  and  CT  are  known,  a„,  U„,  u  and  v  are 
e  e  e  L  Hop 

uniquely  determined.  The  resulting  equations  of  state  for 
aluminum  at  20°C  are  given  in  Figures  13,  14  and  15  and  are 
listed  below: 


=  1.0  +  140.4  u  +  37.7  u  (Std.  Error  0.2  kbar) 
P  P 

=  2.3  +  708  y  +  1607  y2 
=  708  y+  1607  y2  (Hydrostat) 


=  5.24  +  1.40  u 


where  y  =  (v  /v' 


These  equations  of  state  are  based  on  the  measurements  indi¬ 
cated  by  the  symbols  in  Figure  13  and  on  a =  5.7  kbar,  and 
are  considered  to  be  valid  over  the  stress  range  of  8  to  200 
kbar.  Accuracy  of  the  data  points  is  estimated  to  be  ±  2% 

for  cj„  and  ±  1%  for  u  . 

H  p 

The  hydrostat  (P  -y)  was  obtained  by  subtracting  the  stress- 
xi 

axis  intercept  value  from  the  a^-y  equation.  This  approach 
satisfies  the  requirement  that  the  P  -y  function  pass  through 
the  origin  and  assumes  the  deviatoric  stress  (o^)  is  constant 
at  2.3  kbar  (no  strain  hardening).  Alternative  ways  of  esti¬ 
mating  are  discussed  in  Appendix  D. 


The  low-pressure  hugoniot  of  6061-T6  was  reported  previously  by 

(51) 

Lundergan  and  Herrmann  ,  who  used  the  flat-plate  impact  tech¬ 
nique  and  free-surface  velocity  measurements  to  determine  the 
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equation  of  state  up  to  31  kbars .  Results  were  given  as  a 
linear  elastic  portion  (aH  =  1044  (y/l+y) )  up  to  6.4  kbar  and 
a  linear  plastic  portion  from  6.4  to  31  kbar  represented  by: 

cK  =  1.6  +  794  (y/l+y)  (12) 


Comparing  this  with  the  hugoniot  obtained  in  the  present  study 
gives : 


y 

aR  (Eq.  11) 
kbar 

0R  (Eq. 
kbar 

0.01 

9.5 

9.5 

0.02 

17.1 

17.2 

0.03 

25.0 

24.7 

0.04 

33.2 

32.1 

f 

The  agreement  is  good  over  the  stress  range  reported  by 
Lundergan  and  Herrmann. 

* 

The  temperature  dependence  of  the  hugoniot  for  6061  was  deter- 

(52) 

mined  using  x-cut  quartz  gages  .  Measurements  were  made  at 
300°  and  550°C  and  results  are  given  in  Figure  16.  The  stress- 
particle  velocity  hugoniots  are  given  by: 


* 

Although  the  starting  material  was  6061-T6,  the  heating  times 
were  long  enough  such  that  diffusion  processes  were  irrever¬ 
sibly  changing  material  properties  and  the  "-T6"  designation 
is  no  longer  applicable. 
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300°C:  =  0.5  +  127  u  +  49.4  u  2 

H  P  P 

(Six  data  points) 

550°C:  a„  =  113.5  u  +  67  u  2 

H  P  P 

(Five  data  points) 

(a„  <  50  kbar) 

n 

ULTRASONIC  EQUATIONS  OF  STATE 

The  equation  of  state  can  also  be  determined  from  ultrasonics 
measurements  (see,  e.g..  Refs.  53-57).  The  very  accurate  mea¬ 
surement  of  elastic  wave  velocities  (and,  therefore,  bulk  rrr  lu¬ 
lus)  as  a  function  of  hydrostatic  pressure  makes  it  possible 
to  directly  estimate  shock  wave  compression  behavior  or  to  com¬ 
bine  ultrasonics  data  with  shock  wave  data  to  predict  very  high 
pressure  response.  Numerous  analytical  and  empirical  relations 
have  been  developed  for  relating  pressure,  volume  and  bulk  modu¬ 
lus,  including  those  of  Birch^®'5^,  Murnaghan  and 
Keane ^2^ .  The  use  of  these  relations  is  discussed  in 
Appendix  E  and  results  are  summarized  below: 

Birch  EOS 
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Murnaghan  EOS ,  Adiabat 


PS  =  158 


f  4.60  1 

.3  I  (l+y)  -lj 


Keane  EOS,  K  =1.93 


ps  =  896 


1.93  "I 

(1+m)  -lj 


-  1004  In  (.1  fy) 


In  the  above  equations,  subscripts  T  and  S  indicate  isothermal 
and  isentropic,  respectively. 


Duvall  gives  a  correction  for  the  change  in  entropy  across  a 
shock  front  based  on  the  Murnaghan  equation  of  state: 


Y  KS 

p  =  p  +  lo 

FH  *S  12 


PH  ~  PS  +  676 


(4  * l)  (l.  ) 

(*)’ 


Zel'dovich  gives  an  expression  relating  P_  and  P_,  assuming 


y/v  ir  constant: 


(56,64) 


Ps  =  PT  +  6TqK0  exp  ;  y, 


(*))- 


(l+u) 


Ps  =  PT  +  14.2  exp  J  1.97  (  x~  )  ~  (l+u) 
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The  above  equations  of  state  are  compared  to  the  shock  wave 
EOS  in  Table  III.  The  entropy  correction  based  on  the  Murnaghan 
equation  has  been  applied  to  the  Birch  and  Keane  equation ;  as 
well.  The  ultrasonic  equations  of  state  show  good  agreement 
with  each  other  and  with  the  shock  wave  hydrostat,  with  the 
maximum  deviation  being  less  than  ±  2%  up  to  115  kbars . 

The  shock  velocity-particle  velocity  relationship  can  also  b^ 
determined  from  ultrasonic  data.  Following  the  methoi  of 
Ruoff^^  one  obtains: 

U  =  Cn  +  Su  +  Au  2  '19) 

S  B  p  p 

where  C_ 

B 

s 

A 


For  aluminum: 

U®  5.19  +  1.40  u  -  0.12  u2 
S  P  P 

The  u2  coefficient  (A  *  -0.12)  was  calculated  for  K~c  =  -0.02. 
If  one  assumed  =  0/  then  A  =  +  0.12.  Since  the  uncertainty 

S  m  U  D 

in  Kqs  is  ^  ±  0.02,  the  ’A!  coefficient  will  be  taken  as  zero, 
which  gives: 

Uc  =  5.19  +  1.40  u  (21) 

b  p 
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This  is  compared  with  the  linear  Ug  -  relation  obtained 
from  low-pressure  shock  wave  data  in  Figure  15. 

Vaidya  and  Kennedy  measured  the  hydrostatic  compressibility 
of  6061-T6  aluminum  up  to  45  kbar.  Results  were  given  as: 

■  (r  =  =  i-2281  x  l°_3pT  -  8.793  x  10-V 

o  ^  L  (22) 

+  0.10266  x  10_6P^ 

Comparing  this  with  the  shock  wave  and  ultrasonic  results  pre¬ 
sented  in  this  report  gives: 


* 


u 

PT(Eq.  11) 
kbar 

pT(Eq.  14) 
kbar 

pT(Eq. 

kbar 

0.0125 

8.9 

9.0 

10.8 

0.0250 

18.4 

18.4 

22.5 

0.0375 

28.2 

28.2 

34.5 

0.0500 

38.6 

38.5 

45.7 

* 

Corrected  for  isothermal  conditions. 


The  hydrostatic  data  of  Vaidya  and  Kennedy  are  'v  20%  above 
(PT  -  u  plane)  the  shock  wave  and  ultrasonic  data.  The  rea¬ 
son  for  this  large  difference  is  not  known. 

YIELD  BEHAVIOR 

A  complete  description  of  the  low-pressure  equation  of  state 
requires  consideration  of  yield  behavior  or  the  hugoniot  elas¬ 
tic  limit.  Values  of  compressive  yield  were  obtained  by  three 
independent  methods. 
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A.  Uniaxial  Stress 

The  yield  level  was  determined  in  uniaxial  stress 
(o  )  as  a  function  of  strain  rate  and  then  con¬ 
verted  to  uniaxial  strain  using: 


Strain  Rate 
-1 


sec 


10 

10 


3 

-  3 


B.  Wave  Profiles 


0 

e 

kbar 

5.8 

5.7 


The  elastic  precursor  level  was  measured  as  a  function 
of  propagation  distance  using  quartz  gages.  These 
profiles  are  discussed  in  the  section  on  Wave  Pro¬ 
pagation. 


Propagation 

Distance 

mm 


12 


(12.7  mm  Stock) 
'  e 

_ kbar _ 

7.5 

5.5 


C.  Hugoniots 

The  yield  level  in  uniaxial  strain  was  inferred  by 
comparison  of  the  elastic  (. qC^)  and  plastic 
( *,  vs.  u  )  hugoniots,  i.c.,  by  defining  the  clastic 

ii  p 
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limit  as  the  intersection  of  the  elastic  response 
line  and  the  curve-fit  to  the  stress-particle  velocity 
data 

10  paCtu  =  A  +  Bu  +  Cu  2  (=  a  |  (24) 

0  L  e  e  e  \  e) 

Solving  for  ua  gives 

u  =  0.033  mm/jisec 
e 

o  =  5.7  kbar 
e 

For  an  isotropic,  polycrystalline  metal  with  elas tic-perfectly 
plastic  behavior  and  no  strain-rate  or  time-dependent  effects, 
one  would  expect  the  same  yield  to  be  obtained  by  each  method. 
Although  not  strain-rate  sensitive  in  the  uniaxial  stress  tests, 
6061-T6  does  show  relaxation  of  the  elastic  wave  in  uniaxial 
strain.  However,  the  elastic  wave  approaches  an  "equilibrium" 
yield  of  5-6  Kbar  which  is  essentially  the  same  as  the  yield 
obtained  from  uniaxial  stress  or  hugoniot  data.  An  average 
value  of  5.7  kbar  was  used  in  determining  the  au  -  p  and 

n 

U  -  u  equations  given  above  (Equation  11)  for  6061-T6  in 
b  P 

the  low-pressure  region. 
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SECTION  V 

WAVE  PROPAGATION 


Compressive  wave  development  and  elastic  precursor  decay  are 
used  as  basic  input  data  for  material  response  models,  while 
a  complete  understanding  of  shock  wave  propagation  requires 
inclusion  of  release  wave  behavior  (see,  e.g.,  Refs.  23, 

67-70) .  Profiles  of  shock  waves  propagated  through  a  specimen 
are  obtained  as  stress-time  or  veloci ty-time  histories.  Pro¬ 
files  may  be  used  in  this  form  directly  for  comparison  with 
calculated  or  predicted  material  response.  When  it  is  neces¬ 
sary  to  transform  the  data  to  obtain  material  stress  from 
transducer  stress  or  particle  velocity  from  interface  velocity, 
an  impedance  matching  technique  is  applied,  utilizing  the  con¬ 
servation  relations  and  assuming  time-independent  behavior. 

The  implications  of  this  assumption  must  be  considered  when 
utilizing  the  transformed  data.  All  wave  profiles  obtained 
(quartz  gage  and  velocity  interferometer)  are  given  in  Appen¬ 
dix  F  for  reference. 


COMPRESSIVE  WAVE  BEHAVIOR 


Structure  in  the  compressive  wave  in  6061-T6  is  shown  in  the 

* 

quartz  gage  data  in  Figure  17.  The  wave  front  is  character¬ 
ized  by  an  elastic  portion  and  a  transition  to  a  spreading 

* 

In  this  and  subsequent  figures  containing  quartz  gage  data, 
tilt  refers  to  the  time  required  for  a  step-input  to  sweep 
across  the  gage  electrode  diameter. 
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U  -  5.  24  +  1.40  u 
S  P 


MATERIAL:  6061-T6AL,  12.7mm  PLATE 
TEST  NUMBER:  1  2  3 

VELOCITY  (mm/ns):  0.257  0.262  0.272 

THICKNESS  (mm  ):  2.01  6.06  11.77 

TILT  (ns):  9  21  19 


.on  .20  .40  .60  .10  1.00  1 

TIME  IN  MICRO  SECONO 


Figure  17  Compressive  Waves,  Propagation 
Distance  Dependence 
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plastic  wave  (nominal  final  stress  of  20  kbar) .  (The  elastic 
wave  would  be  overdriven  by  the  plastic  wave  at  ^  140  kbar.) 

The  elastic  portion  is  not  appreciably  influenced  by  peak 
stress,  at  least  for  stresses  up  to  4  times  the  elastic  limit, 
as  shown  in  Figure  18.  Also  shown  in  Figures  17  and  18  are 
shock  velocities  as  calculated  from  the  shock  wave  hugoniot 
(Equation  11).  Although  there  is  qualitative  agreement  in 
that  the  predicted  plastic  wave  velocity  lies  within  the  rise¬ 
time  of  the  measured  wave,  the  plastic  wave  cannot  be  considered 
to  have  a  unique  velocity  for  a  given  stress  level,  at  least 
for  stresses  below  s,  25  kbar.  This  means  that  transformation 
of  the  hugoniot  from  the  o„  -  u  plane  to  ou  -  u  or  U„  -  u 

H  p  H  bp 

by  assuming  ideal  elastic-plastic  wave  structure  and  steady- 
state  conditions  may  lead  to  some  uncertainty  at  lower  stresses. 

Stress-particle  velocity  paths  were  calculated  from  quartz  gage 
transmitted  wave  profiles  and  are  compared  to  the  hugoniot  in 
Figure  19.  The  calculations  were  made  using  QZ,  a  program  which 
assumes  steady-state  behavior  and  applies  the  hugoniot  jump 
conditions  incrementally  to  the  wave  profile  by  an  impedance 
matching  technique . ^ 50 ^  The  dynamic  paths  generally  lie  above 
the  equilibrium  hugoniot  with  only  the  end  points  or  final  states 
approaching  the  equilibrium  values.  This  deviation  is  due  to 
time-dependent  effects. 

The  conditions  for  development  of  steady-state  plastic  waves  in 
6061-T6  were  studied  by  conducting  a  series  of  velocity  inter¬ 
ferometer  tests  and  combining  the  results  with  similai  data 
reported  by  Johnson  and  Barker.  J  In  an  clastic-plastic  wave 
system,  the  effects  of  time-dependency  and  non-linearity  in¬ 
itially  car.  cause  the  plastic  portion  to  spread  as  it  propa¬ 
gates.  Eventually,  however,  the  plastic  portion  may  achieve 
a  stable  condition  and  propagate  relatively  unchanged  (provided 
conditions  of  u.iiaxial  strain  arc  maintained  and  the  wave  front 
is  not  overtaken  by  a  release  wave) . 
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Compressive  wave  profiles  are  shown  in  Figure  20  for  two  series 
of  tests  using  the  velocity  interferometer:  varying  peak  stress 
at  constant  propagation  distance,  and  varying  propagation  dis¬ 
tance  at  constant  peak  stress.  (Test  conditions  for  velocity 
interferometer  data  discussed  in  this  report  are  listed  in 
Table  IV.  The  test  numbers  correspond  to  those  listed  next  to 
each  wave  profile  in  subsequent  figures.)  At  26  kbar  Jdig 
rise-time  in  the  plastic  wave  increases  between  3  and  12  mm, 
but  wave  stability  is  achieved  by  20  mm.  The  series  at  3  mm 
thickness  shows  a  steepening  plastic  wave  with  increasing  peak 

stress.  A  more  complete  mapping  of  the  formation  of  steady- 

* 

state  waves  is  given  in  Figure  21.  Tests  at  a  given  stress 
which  showed  steady-state  behavior  in  the  plastic  wave  are 
shown  by  solid  symbols.  The  shaded  area  shows  the  indistinct 
border  between  non-steady-state  and  steady-state  behavior. 

Dislocation  models  for  predicting  compressive  wave  development 
utilize  measurements  of  elastic  precursor  decay  to  evaluate 
various  constants  related  to  dislocation  density  and  multipli¬ 
cation.  72^  Unfortunately,  it  is  very  difficult  to  obtain 
accurate  elastic  precursor  level  data  at  propagation  distances 
less  than  'v  2  mm  and  it  is  in  this  region  (0  to  2  mm)  that  most 
of  the  precursor  decay  occurs  for  many  metals  of  interest. 

Also,  even  at  larger  distances  there  is  always  some  uncertainty 
in  the  precursor  level  due  to  both  material  scatter  and  disper¬ 
sion  characteristics  and  to  experimental  errors  (particularly 
tilt) .  The  result  is  that  much  of  the  data  presently  available, 
when  given  realistic  uncertainty  bars,  can  be  used  to  support 
a  number  of  different  models  of  precursor  attenuation.  The 

* 

Steady-state  conditions  were  studied  by  comparing  rise-times 
in  the  plastic  wave  at  approximately  20  to  80%  stress  levels 
above  the  elastic  limit. 
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precursor  decay  data  for  6061-T6  is  summarized  in  Figure  22. 

The  calculated  initial  elastic  impact  stress  was  23  kbar  for 
these  tests  and  the  precursor  level  has  dropped  ^  70%  in 
2  mm  of  travel.  Over  the  range  of  propagation  distance  studied, 
the  precursor  level  in  12.7  mm  stock  was  ^  1  kbar  greater  than 
in  the  3.2  or  6.4  mm  stock. 


Figure  20  Compressive  Wave  Development,  6061-T6 
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TABLE  IV 

VELOCITY  INTERFEROMETER  TEST  DATA,  6061-T6  ALUMINUM 


VI 

o 

X0 

X 

X/XQ 

Teat 

No. 

Impact 

Velocity 

(ram/us) 

Max. 

Stress 

(kbar) 

Impactor 

Thickness 

(mm) 

Target 

Th^ckpess 

Stock 

(mm) 

IS 

0.269 

21 

1.879 

3.124 

1.66 

3.2 

50 

0.222 

17 

1. 552 

3.081 

1.99 

3.2 

52 

0.2  52 

19 

1.587 

12.220 

7.70 

12.7 

53 

0.284 

22 

1.440 

18.192 

12.63 

50.8 

612 

0.510 

39 

1.580 

3.106 

1.97 

3.2 

63 

0.193 

15 

1.570 

3.058 

1.95 

3  2 

732 

1.242 

103 

0.269 

2.520 

9.37 

3.2 

742 

1.238 

102 

0.252 

5.067 

20.15 

6.4 

792 

1.235 

102 

0.252 

1.029 

4.09 

1.6 

83 

0.261 

20 

0.475 

1.181 

2.49 

3.2 

104 

0,255 

19 

.590 

3.080 

1.94 

3.2 

108 

0.263 

20 

1.562 

2.776 

1.78 

3.2 

109 

0.340 

26 

0.431 

3.985 

9.23 

6.4 

110 

0.418 

32 

0.419 

3.970 

9.47 

6.4 

133 

0.265 

20 

1.600 

3.073 

1.92 

3.2 

135 

0.247 

19 

1.598 

3.17B3 

1.99 

1.6 

143 

0.306 

23 

0.904 

3.061 

3.39 

3.2 

177 

0.275 

19 

2.982 

3.249 

1.09 

50.8 

178 

0.373 

25 

3.058 

2.972 

0.97 

50.8 

181 

0.358 

24 

1.0!8 

6.033 

2.00 

50.8 

182 

0.354 

24 

3.041 

12.141 

3.97 

50.8 

183 

0.329 

25 

12.649 

20.259 

1.60 

50.8 

187 

0.489 

43 

3.157 

2.977 

0.94 

50.8 

1.  a.g.,  Al-Al  is  aluaainuni  iapactor  into  aluminum  target 
with  Iraa  raar  air  face  (FO  ia  fused  t,u»rii  window). 

2.  Plaatl^laa  -  backed  inpact or. 

).  Artificial  spall  plana  at  X/2. 


Conflg.1 


Al-Al 

Al-Al 

Al-Al 

Al-Al 

A1-A1/FQ 

Al-Al 

A1-A1/FQ 

A1-A1/FQ 

Al-Al/FQ 

Al-Al 

Al-Al 

Al-Al 

Al-Al 

Al-Al 

Al-Al 

Al-Al 

Al-Al 

ro-Ai/ro 

ro-Ai/ro 

ro-Ai/ro 

ro-Ai/ro 

Al-Al/TO 

ro-Ai/ro 
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Figure  21  Steady  State  Wave  Conditions,  6061-T6 


the  uncertainty  bars  in  Figure  22  ref.lect  the  dispersion  and 
rounding  at  the  elastic  front.  This  spreading  is  due  to  at 
least  two  factors.  First,  since  the  quartz  gage  averages  stress 
over  the  electrode  area,  small  differences  in  wave  front  arrival 
times  at  the  specimen/gage  interface  would  give  a  finite  rise¬ 
time  rather  than  an  instantaneous  stress  jump.  These  snail 
differences  in  arrival  times  may  be  due  to  wave  reflections  and 
dispersion  at  grain  boundaries  and  to  slight  differences  in 
local  wave  velocity  arising  from  grain  anisotropy.  Second,  the 
influence  of  shock  wave  tilt  on  a  finite-area  gage  is  to  smooth 
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out  abrupt  changes  in  stress  level  as  well  as  to  increase  re¬ 
cording  rise- time  in  the  wave  front.  This  becomes  particularly 
severe  at  short  propagation  distances  where  separation  of  the 
elastic  and  plastic  waves  is  only  on  the  order  of  tens  of  nano¬ 
seconds  . 
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Figure  22  Elastic  Precursor  Decay,  6061-T6 
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Of  about  15  quartz  gage  transmitted  wave  tests  conducted  with 
6061-T6,  only  one  showed  relaxation  behind  the  elastic  precursor, 
as  shown  in  Figures  23  and  24.  This  relaxation  was  observed  for 
a  propagation  distance  of  11.25  mm  with  50.8  mm  stock,  but  was 
not  observed  for  propagation  distances  of  6  or  20  mm  in  targets 
machined  from  the  same  plate.  Tilt  was  slightly  higher  for  the 
6  mm  test  (23  nsec) ,  which  may  have  smeared  the  wave  front  such 
that  relaxation  could  not  be  resolved.  For  the  20  mm  test,  which 
had  a  tilt  of  16  nsec,  the  precursor  "spike"  may  have  relaxed 
completely.  Material  variability  could  also  have  influenced 
the  results. 


12.7  mm  PLATE  50.8  mm  PLATE 


AI-»AI/Q  AI-*AI/Q 

Vj«  0.272  rnm/ps  V.-  0.271  mm/ps 

CTs-  20.8  kbar  <3-  20.7  kbar 

Figure  23  Quartz  Gage  Transmitted  Wave 

Records,  6061-T6 
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Figure  24  Compressive  Waves,  Material  Differences 
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Compressive  wave  structure  was  also  measured  in  annealed 
6061,  and  is  compared  with  6061-T6  in  Figure  25.  The  elastic 
limit  has  decreased  to  ^  1  kbar  for  this  propagation  distance 
(3.1  mm) . 


§ 


Figure  25 


Compressive  Waves,  Annealing  Effect 
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RELEASE  WAVES  AND  WAVE  ATTENUATION 

Release  waves  in  6061-T6  were  studied  by  using  relatively  thin 
impactors  and  measuring  the  complete  wave  profile  with  the 
laser  velocity  interferometer.  After  impact,  elastic  and 
plastic  waves  propagate  from  the  impact  surface  into  both  im- 
pactor  and  target.  The  waves  that  reach  the  impactor  rear 
surface  reflect  as  rarefaction  waves,  unloading  the  material. 

These  unloading  waves  propagate  across  the  impact  surface  into 
the  target  and  are  recorded  at  the  target  rear  surface  after 
arrival  of  the  compressive  wave. 

Three  unattenuated  wave  profiles  for  6061-T6  are  shown  in 
Figure  26,  one  with  a  6061-T6  impactor  and  two  with  fused 
quartz  impactors.  (Additional  impact  data  are  given  in  Table  IV) . 
All  three  targets  had  fused  quartz  windows  on  the  rear  face  to 
reduce  the  magnitude  of  unloading  due  to  reflection  of  the  com¬ 
pressive  wave.  This  prevents  spall  fracture  in  the  target. 

The  6061-T6  impactor  for  test  61  was  backed  by  Plexiglas  to 
prevent  bowing  during  launch  and  this  accounts  for  the  wave  not 
releasing  completely.  Also,  the  measured  interface  velocity  for 
these  tests  is  slightly  higher  than  particle  velocity  in  the  tar¬ 
get  since  the  impedance  of  6061-T6  is  slightly  greater  than  fused 
quartz.  Structure  is  evident  in  the  release  wave  for  all  three 
tests,  however,  the  development  of  this  structure  for  test  61 
with  a  6061-T6  impactor  is  more  complex  than  for  tests  with  fused 
quartz  impactors.  This  is  because  elastic  and  plastic  waves 
spread  in  the  6061-T6  impactor  such  that  all  the  unloading  does 
not  occur  at  the  same  time,  i.e.,  the  release  wave  does  not 
start  out  as  an  abrupt  pressure  drop.  The  fused  quartz  impactor 
gives  a  "square"  input  wave  to  the  target,  with  a  nearly  instan¬ 
taneous  release  at  the  impactor- target  interface.  Therefore, 
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the  structure  in  the  release  waves  for  tests  181  and  182  should 
be  due  to  elastic-plastic  unloading  in  the  target  material  only. 
A  shock  wave  attenuation  study  by  McMillan,  Isbell  and  Jones 
indicated  that  the  initial  pressure  release  in  6061-T6  is  elas¬ 
tic  even  at  pressures  up  to  1  mbar  and  the  elastic  release  wave 

(73) 

amplitude  increases  with  increasing  peak  pressure. 

If  target  thickness  is  great  enough  compared  to  impactor  thick¬ 
ness,  the  release  wave  will  overtake  the  compressive  wa\  and 
attenuate  the  peak  stress.  Wave  attenuation  is  shown  in  Fig¬ 
ure  27  for  an  initial  stress  of  20  kbar  and  in  Figure  28  for 
an  initial  stress  of  102  kbar.  At  the  higher  stress  level,  the 
impactors  were  backed  by  Plexiglas  so  the  stress  wave  was  not 
released  to  zero.  The  attenuation  results  are  shown  in  Fig¬ 
ure  29  in  terms  of  decrease  in  peak  stress  vs.  X/X^ .  Also  shown 
are  calculated  overtaking  points  assuming  ideal  elastic-plastic 
compression  and  release  waves.  This  is  given  by: 


(25) 


release  wave  velocity  (6.5  mm/us 
mm/us  for  102  kbar^7^). 

for  102  kbar  are  consistent  with 
the  calculated  overtaking  point,  but  the  results  for  20  kbar 
show  attenuation  much  sooner  than  calculated.  This  can  be 
attributed  to  the  ramped  plastic  compressive  wave  at  this 
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where  C?-  is  elastic 

ii 

for  20  kbar  and  7.0 


The  experimental  results 
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SECTION  VI 

SPALL  FRACTURE 


Spall  fracture  by  plate  impact  results  from  reflection  of 
compressive  waves  from  a  relatively  low  impedance  interface 
(normally  a  free  surface)  and  subsequent  wave  interaction 
(see,  e.g..  Refs.  74-76).  Spall  studies  were  carried  out 
with  both  passive  and  active  techniques.  Passive  methods 
involve  the  recovery  and  examination  of  shock- loaded  speci¬ 
mens.  Metal lographic  examination  establishes  the  type  and 
degree  of  damage,  which  can  bs  correlated  with  test  parameters 
such  as  velocity,  impactor  thickness  and  target  thickness. 
Active  methods  utilize  the  laser  velocity  interferometer  and 
provide  quantitative,  time-resolved  data  on  the  influence  of 
internal  fractures  or  spall  surfaces  on  shock  wave  profiles, 
as  measured  at  the  rear  surface. 


RECOVERY  TESTS 


The  spall  behavior  of  6061-T6  was  studied  by  carrying  out  a 
series  of  impact  and  recovery  tests,  where  the  target  was 
sectioned  across  a  diameter,  polished,  etched,  and  examined 
optically  at  a  magnification  of  50X.  The  specimens  were  then 
graded  or  classified  according  to  the  degree  of  fracture  that 
was  observed,  which  ranged  from  no  visible  damage  to  complete 
material  separation.  The  onset  of  significant  fracture  is 
generally  referred  to  as  incipient  spall  and  is  a  critical 
point  since  it  can  be  used  to  deduce  the  dynamic  fracture 
strength,  and  corresponds  to  generation  of  sufficient  free 
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surfaces  within  the  material  to  reflect  a  portion  of  the 
interacting  release  waves  as  a  compressive  wave.  The  incip¬ 
ient  spall  threshold  was  defined  as  the  impact  velocity  (for 
a  given  set  of  impact  parameters)  corresponding  to  cracking 
over  at  least  50%  of  the  width  of  the  section  estimated  to  be 
under  a  condition  of  plane  strain  during  the  time  of  loading. 


Incipient  spall  results  from  all  tests,  including  an  attenuated 
pulse  series  and  elevated  temperature  tests,  are  summarized  in 
Table  V.  The  data  for  all  20°C  tests  with  unattenuated  pulses 
(X/Xq  =»  2)  are  given  in  Figure  30.  The  symbols  are  defined  as: 

•  Complete  Separation 
0  Above  Incipient 

•  Incipient  Spall 
e  Below  Incipient 

0  No  Visible  Damage 


As  is  typically  the  case  with  metals,  the  impact  velocity  for 

incipient  spall  increases  for  decreasing  impactor  thickness. 

This  implies  that  peak  stress  (and,  therefore,  strain)  required 

to  create  spall  fractures  increases  with  decreasing  pulse  width 

or  time  of  loading.  These  results  are  compared  to  those  of 

Kreer'  and  Butcher,  et  al.  'in  Figure  31.  These  data  are 

* 

all  for  gun-launched,  plate-impact  techniques.  Although 
Butcher's  tests  were  for  X/XQ  *  4,  a  comparison  with  data  for 

l - - - 

Spall  in  6061-T6  has  been  studied  by  Keller  and  Tuler(78), 
Jajosky  and  Ferdman(79)  and  Cohen  and  Berkowitz (80) f  using 
exploding-foil  techniques;  however,  their  results  are  not 
included  here  since  differences  in  experimental  method,  con¬ 
dition  of  the  impactors  and  metallographic  analysis  prevent 
a  direct  comparison.  Also,  Stefansky  and  Shea (81)  studied 
dynamic  fracture  in  several  metals,  including  6061-T6,  using 
electron  beam  pulses. 
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X/XQ  =  2  is  valid  since  the  incipient  spall  threshold  in 
6061-T6  increases  only  ^  5%  for  a  change  in  X/Xq  =  2  to 
X/Xq  =  4,  as  shown  in  Figure  32.  The  data  given  in  Figure  31 
show  reasonable  agreement  for  XQ  >  0.8  mm.  For  thinner  impactors , 
however,  the  data  from  this  work  and  that  of  Kreer  deviate  signi¬ 
ficantly.  This  may  be  due  to  differences  in  experimental  tech¬ 
niques  (e.g.,  impactor  bowing,  tilt,  target  construction  or 
metallographic  methods) ,  although  Kreer  indicates  that  material 
property  variation  may  also  be  a  cause. 

Optical  photomicrographs  of  recovered  specimens  are  shown  in 

★ 

Figures  33  to  35.  The  specimens  were  polished  and  etched  to 
provide  maximum  contrast  between  cracks  and  sound  material. 

The  procedure  followed  was : 

1.  Wet  grind,  320/400/600  grit  silicon  carbide  paper. 

2.  Rough  polish,  5  ym  aluminum  oxide  slurry  on 
rayon  cloth. 

3.  Final  polish,  0.05  ym  gamma  alumina  slurry  on 
rayon  cloth. 

4.  Etch,  10  ml  hydrofluoric  acid  (HF) ,  15  ml  hydro¬ 
chloric  acid  (HCl)  and  90  ml  water,  30-90  sec 
immersion  (20°C),  wash'in  warm  water. 

As  noted  above,  spall  cannot  be  rigorously  defined  in  terms  of 
a  unique  impact  velocity  or  stress,  but  requires  consideration 
of  the  degree  of  fracture.  In  general,  as  impact  velocity  in¬ 
creases,  spall  progresses  from  isolated  cracks  or  voids  only 
a  few  grains  in  size  to  complete  material  separation  across  the 
entire  spa3 1  plane. 

* 

All  spall  photomicrographs  in  this  report  are  oriented  such 
that  initial  shock  wave  propagation  was  from  bottom  to  top. 
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Fracture  in  6061-T6  is  of  a  ductile  nature  with  the  development 
of  individual  voids.  As  velocity  increases  (see,  e.g.  Fig¬ 
ures  33  and  34)  the  voids  begin  to  coalesce  and  a  complete 
fracture  surface  is  formed.  Note  that  there  is  some  disper¬ 
sion  of  the  voids  around  a  nominal  "spall  plane".  This  was 

( 82 )  /  o  3 \ 

also  found  for  titanium  and  tantalum'  ,  which  have  well- 
defined  elastic  wave  structures,  while  copper has  a  poorly- 
defined,  low-amplitude  elastic  wave  and  showed  very  little  void 
dispersion  around  the  spall  plane.  Spall  fractures  in  the  four 
metals  are  compared  in  Figure  36 . 

Scanning  microfractographs  of  the  fracture  surfaces  for  completely 
separated  specimens  are  shown  in  Figures  37  and  38.  The  voids  ap¬ 
pear  to  have  initiated  at  brittle  particles  (precipitates  or 
second-phase  particles)  which  are  visible  in  the  center  of  the 

voids.  Some  of  the  particles  themselves  appear  to  be  fractured, 

( 85 ) 

an  observation  previously  report  by  Au  for  6061-T6.  Complete 
fracture  has  been  primarily  by  normal  rupture  with  evidence  of 
shearing  near  the  necked  or  peripheral  areas  of  the  voids.  There 
was  essentially  no  difference  in  fracture  surface  for  specimens 
taken  from  different  stock  (see  Figure  38).  Tuler'  reported 
essentially  no  change  in  fracture  surface  in  6061-T6  for  failure 
due  to  slow  bending,  charpy  impact  and  plate  impact. 

Several  tests  were  conducted  to  determine  the  effect  of  preshock¬ 
ing  on  spall  behavior.  A  target  was  assembled  with  radial  and 
rear  momentum  traps  and  subjected  to  a  20  kbar,  0.5  nsec  compres¬ 
sive  pulse.  The  recovered  specimen  was  then  tested  at  a  velo¬ 
city  which  caused  spall  in  a  specimen  that  had  not  been  pre¬ 
shocked.  Results  are  compared  in  Figure  39.  Apparently,  pre¬ 
shocking  has  increased  the  spall  fracture  strength,  since  the  un¬ 
shocked  specimen  shows  spall  fractures  while  the  preshocked  spec¬ 
imen  does  not,  for  the  same  test  conditions. 
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NOT  PRESHOCKED  PRESHOCKED  TO  20  kbars 


Figure  39  Spall  Fractures,  Preshocking  Effect, 
6061-T6  (20°C) 


Stevens  and  Tuler^87^  studied  the  influence  of  varying  amounts 
of  precompression  on  spall  fracture  by  backing  the  target  with 
a  material  of  lower  impedance  (rather  than  a  free  rear  surface) . 
They  reported  that  increasing  the  precompression  stress  from 
12  to  82  kbar  did  not  significantly  change  the  calculated  in¬ 
cipient  spall  tensile  stress. 

Temperature  dependence  of  spall  for  6061  is  given  in  Figure  40. 
The  impact  velocity  for  incipient  spall  decreases  with  increas¬ 
ing  temperature.  This  behavior  is  probably  related  to  changes 
in  the  dynamic  fracture  strength  of  the  material.  However,  as 


77 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 


MSL-70-23,  Vol.  Ill 

Charest  has  pointed  out^88\  the  use  of  a  room  temperature  im- 
pactor  and  a  hot  target  (as  was  done  here) ,  gives  impedance 
differences  that  influence  development  of  the  release  wave 
structure  and  therefore  stress-time  history  in  the  target. 

Also,  the  reduced  hugoniot  elastic  limit  in  6061  at  elevated 
temperature  will  affect  the  wave  structure  in  the  target. 

Optical  photomicrographs  and  scanning  microfractographs  of  6061 
spall  specimens  tested  at  several  temperatures  are  shown  in 
Figures  41  and  42.  Although  the  impact  velocity  required  for 
incipient  spall  and  for  complete  separation  has  decreased, 
there  has  been  no  appreciable  change  in  the  appearance  of  voids 
or  fracture  surface. 

Several  specimens  were  annealed  by  heating  at  400°C  for  three 
hours  and  then  tested  at  20 °C.  The  incipient  spall  threshold 
was  ^  0.145  mm/ysec  or  ^  25%  less  than  6061-T6  for  the  same 
thickness.  This  is  comparable  to  the  difference  reported  by 

(CO) 

Butcher  for  "-T6"  and  "-0"  condition  material.  Optical 
photomicrographs  of  6061-0  spall  specimens  are  shown  in  Fig¬ 
ure  43,  and  spall  fracture  surfaces  of  6061-0  and  6061-T6  are 
compared  in  Figure  44. 


SPALL  PROFILES 

Spall  wave  profiles  were  studied  using  an  aluminum-into-aluminum 
configuration,  with  a  free  rear  surface  on  the  target  for  com¬ 
plete  reflection  of  the  compressive  front.  A  schematic  of  the 
wave  interactions  resulting  from  such  a  test  in  shown  in  Fig¬ 
ure  45.  Several  simplifying  assumptions  are  made.  The  shock 
is  shown  as  a  single  wave  and  the  release  waves  originating  at 
boundaries  are  considered  simple,  centered  rarefaction  fans. 
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Figure  43  Spall  Fractures,  6061- 
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WIVE  IITEMCTIOIS  II  lit  VELOCITY  .  TIME  HI5I0IY 

TIME  0IS1IICE  Hill  IT  llil  SUIF1CE 


Figure  45  Wave  Interaction  Schematic,  Spall  Test 

Also,  small  interactions  between  shock  waves  have  been  omitted. 
Fracture  of  the  material  is  assumed  to  occur  along  a  single 
surface  (spall  plane)  rather  than  in  a  volume  around  this  plane 
(see  Figure  33) .  Wave  interactions  at  the  spall  plane  are  shown 
as  if  the  material  exhibits  a  time-dependent  fracture,  i.e., 
the  percentage  of  fractured  area  in  the  spall  plane  increases 
with  time  until  the  entire  spall  surface  has  separated.  During 
the  early  stages  of  the  fracture  process,  before  complete  sep¬ 
aration  is  achieved,  waves  crossing  the  spall  plane  are  part¬ 
ially  transmitted  and  partially  reflected.  After  complete 
separation  has  occurred,  any  waves  trapped  in  the  spalled  piece 
reverberate,  giving  a  sinusoidal  motion  to  the  rear  surface, 
which  "rventually  damps  out  due  to  attenuation  and  dispersion 
in  the  material  and  to  the  effects  of  edge  rarefactions. 
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Measured  profiles  for  unattenuated  pulses  in  6061-T6  are  shown 
in  Figure  46,  with  photomicrographs  of  the  recovered  and  sec¬ 
tioned  targets.  Also  shown  for  reference  are  the  incipient 
and  complete  spall  levels  for  this  impactor/target  combination 
as  determined  from  the  recovery  tests.  Evidence  of  fracture  in 
the  target  is  shown  as  a  reversal  of  the  measured  release  wave. 
The  amount  of  pullback  (decrease  in  free  surface  velocity  to 
the  point  of  first  reversal)  may  be  related  to  the  tensile 
strength  at  the  spall  plane  of  the  material  for  a  given  impact 
geometry.  An  empirical  relation  for  this  is  given  by: 

os  =  pC(Aufs/2)  (26) 

where  o  is  tensile  strength,  p  and  C  are  local 
density  and  longitudinal  wave  speed,  and  Aufg 
is  pullback. 

The  data  in  Figure  46  give  a  *  12  kbar.  This  calculated  spall 
strength  is  not  constant  for  a  g~  ven  material  since  it  is  pro¬ 
portional  to  Aufs  which  has  been  found  to  change  with  such 
factors  as  pulse  width  (see  Figure  47) ,  target  temperature 
(see  Figure  48) ,  and  pulse  shape  (see  Figure  49) .  These  fac¬ 
tors  are  discussed  in  detail  in  Reference  92.  However,  for  a 
given  set  of  impact  conditions  (thickness,  material,  tempera¬ 
ture)  ,  the  pullback  is  relatively  insensitive  to  the  maximum 
compressive  stress  as  well  as  to  the  degree  of  fracture,  as 
shown  in  Figure  50,  This  is  consistent  with  the  varying-pre- 
compression  results  of  Stevens  and  Tuler^®7^  referred  to  ear¬ 
lier. 


Brammer^D  reported  spall  profiles  for  2024  Aluminum  at 
22°C  and  371°C  which  also  showed  a  decrease  in  pullback 
at  the  higher  temperature. 
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In  the  elastic  range,  pullback  equals  impact  velocity,  while 
at  impact  velocities  causing  fracture,  the  pullback  is  approxi¬ 
mately  0.135  mm/ysec  (independent  of  degree  of  fracture  and 
maximum  stress).  Thus,  the  pullback  appears  to  be  associated 
with  the  incipient  spall  threshold,  i.e.,  the  formation  and 
initial  coalescence  of  voids  in  the  material,  rather  than  with 
complete  fracture  or  separation.  This  conclusion  may  allow 
comparison  of  the  dynamic  fracture  strength  of  the  material 
over  a  range  of  wave  shapes  and  times  of  loading  without  the 
necessity  of  conducting  all  tests  with  exactly  the  same  degree 
of  final  fracture.  An  additional  benefit  may  result  if  the 
pullback  can  be  related  to  the  results  of  the  recovery-type 
measurements  of  the  incipient  spall  threshold.  Then  a  single 
wave  profile  test  could,  for  a  given  set  of  initial  conditions, 
provide  the  same  incipient  threshold  data  as  a  series  of  re¬ 
covery  tests.  This  data  could  be  of  use  in  evaluating  various 
spall  fracture  models  through  comparison  of  measured  shock 

( q  3—99 ) 

wave/spall  plane  interactions  with  calculated  wave  profiles. 

The  influence  of  spall  plane  strength  was  studied  with  a  tar¬ 
get  having  an  "artifical"  spall  plane  of  essentially  zero 
strength.  This  zero-strength  plane  was  actually  the  interface 
between  two  target  plates  bonded  together  with  a  thin  (<  5  ym) 
epoxy  layer  having  a  tensile  strength  of  ^  0.01  kbars.  The 
rear  carget  plate  was  the  same  thickness  as  the  impactor, 
putting  the  interface  at  the  nominal  spall  plane.  Referring 
to  Figure  51,  for  a  single  wave  system  the  free  surface  velo¬ 
city  would  be  expected  to  increase  abruptly  to  the  impact  velo¬ 
city.  The  momentum  of  the  impact  is  simply  transferred  to  the 
spall  piece.  For  a  two-wave  system,  the  momentum  transfer  will 
again  take  place.  However,  the  interactions  at  the  interface 
are  now  more  complicated  and  a  net  tensile  stress  may  develop 
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Figure  51  "Zero- Strength"  Spall  Behavior  Schematic 

in  the  spall  piece.  The  pullback  at  the  rear  surface  shows 
evidence  of  this  tensile  pulse.  Figure  52  shows  the  measured 
wave  profile  for  such  a  test.  The  ringing  in  the  spall  piece 
continues  for  many  microseconds.  The  damped  pullback  apparently 
will  decrease  to  zero  amplitude  and  the  piece  will  continue  with 
approximately  impact  velocity. 

Spall  profiles  for  attenuated  pulses,  where  the  time  duration 
of  the  peak  tensile  stress  is  expected  to  be  less  than  for  the 
unattenuated  case,  are  shown  in  Figure  53.  The  pullback  for 
this  case  is  0.175  mm/ysec,  giving  a  spall  strength  of  'v  15  kbar. 


k. 
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Figure  53  Attenuated  Spall  Profiles,  6061-T6 
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SUMMARY 


The  dynamic  properties  of  6061-T6  aluminum  were  measured  under 
uniaxial  stress  and  uniaxial  strain  conditions.  The  material 
tested  was  sheet  and  plate  stock  with  average  hardness  of  57Rn. 

D 

Material  response  to  compressive  uniaxial  stress  deformation 
showed  a  well-defined  yield  with  slight  strain-hardening. 
Strain-rate  sensitivity  was  essentially  zero  at  20°C  (2.9  kbar 
yeild)  but  sensitivty  increased  and  yield  decreased  at  150, 

204  and  288°C,  for  a  range  in  strain  rate  of  0.005/sec  to 
1000/sec. 

Longitudinal  and  shear  wave  velocities  were  measured  as  func¬ 
tions  of  temperature  and  pressure,  giving  the  following  (in 
mm/usec,  ; 


CL  =  6.387  - 
CL  =  3.214  - 
(T  in  °C) 

CT  =  6.368  + 

Xj 

Cs  =  3.197  + 
(P  in  kbar) 


(0.00093T) 

(0.00084T) 

(0.0180P)  - 
(0.0108P)  - 


(0.0001P2) 

(0.00005P2) 


Principal  elastic  constants  were  evaluated  at  r.’  =  20°C  and 
P  *  0,  including: 


Bulk  Modulus  =  Kq  =  728  kbar 
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Shear  Modulus  =  Gq  =  276  kbar 

(#)s  =  2-12'  (#)t  =  2-25 

Poisson's  Ratio  =  0.331 

Elastic  wave  velocity  data  were  also  used  in  determining  the 
Gruneisen  parameter  y  and  the  Debye  temperature  0 : 

y  =  1.97  @  T  =  20°C  and  P  =  0 
|X  =  +  0. 00005 1/°C  a  T  =  20°C 

li-  =  -  0. 0051/kbar  @  P  =  0 
dP 

0  =  420°K 

The  shock  wave  hugoniot  was  determined  experimentally  to  be: 
a  =  1.0  +  140.4  u+  37.7  u2  (a  <  200  kbar) 

tl  p  p  fl 

Assumption  of  an  ideal  elastic-plastic  wave  structure  with  a 
5.7  kbar  hugoniot  elastic  limit  and  equilibrium  initial  and 
final  states  gave: 

ou  =  2.3  +  708  v  +  1607u2 
n 

Pu  *  708  u  +  1607u2 

n 

U_  «  5.24  +  1. 40u 
S  p 

The  hydrostat  as  determined  from  elastic  constants  data  was: 
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PH  =  896 


(l+u) 


1.93 


-1 


-  1004  In  (l+p)  +  676 


(*) 


(Keane  equation  of  state) 


Measurements  of  compressive  wave  development  showed  the  wave 
front  to  be  characterized  by  a  well-defined  elastic  portion 
and  a  transition  to  a  spreading  plastic  wave.  Elastic  pre¬ 
cursor  decay  was  ■v  75%  to  6  kbar  after  6  mm  of  travel.  Peak 
pressure  attenuation  was  ^  44%  at  X/XQ  =  12.6  for  a  20  kbar, 
0.5  psec  pulse,  and  ^  46%  at  X/Xq  =  20  for  a  102  kbar, 

0.08  psec  pulse. 


Spall  fracture  behavior  was  determined  for  several  elastic 
pulse  widths  ranging  from  0.08  to  1.26  psec.  The  incipient 
spall  threshold,  expressed  as  impact  velocity,  increased  with 
decreasing  pulse  width,  and  decreased  with  increasing  tempera¬ 
ture  for  20,  288  and  400°C.  The  spall  plane  was  reasonably 
well-defined,  and  was  located  approximately  one  impactcr  thick¬ 
ness  from  the  target  rear  surface,  for  unattenuated  pulses. 
Failure  was  of  a  ductile  nature  with  voids  initiated  at  rel¬ 
atively  brittle  particles. 
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APPENDIX  A 

ELASTIC  CONSTANTS  EQUATIONS 


Following  are  the  equations  used  in  calculating  various  elastic 

constants  and  the  adiabatic  and  isothermal  pressure  and  tempera- 

* 

ture  derivatives  of  the  bulk  and  shear  moduli i  (K  and  G)  : 

CB  =  V  CL2-4/3  CS2 

(Bulk  Wave  Velocity) 

i - r 

CE  =  y  2  (i+v)  Cg 

(Sound  wave  velocity,  longitudinal 
wave  of  infinite  wave  length  in 
bar,  lateral  displacement  allowed) 

CR:  kj  -  8kJ  4-  (24~16a12)k12  +  (16a12-16)  =  0 

where  k2  =  (CR/Cg)2,  c^2  =  (Cg/CL) 2 
(Rayleigh  or  Surface  Wave  Velocity) 


The  expressions  for  the  pressure  derivatives  are  simplified 
in  these  equations,  e.g.. 
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v  =  [0.5-(Cs/Cl)2]/[1.0-(Cs/Cl)2] 
(Poisson's  Ratio) 


K  =  pC2  =  p(C2-4/3C2) 


(Bulk  Modulus) 


G  =  pc; 


(Shear  or  Rigidity  Modi ius) 


E  =  2p(l+v)  Cg 

(Young's  or  Elastic  Modulus) 


X  =  2p\JCg/(l-2v) 

(Lame's  Parameter) 


T  0 

Ko  "  T  +  6yqt 


where  yn  = 


pocp 


c  <  c  • 

Kos  =  kot 


;•(*)•(¥) 
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APPENDIX  B 

GRUNEISEN  PARAMETER  ESTIMATES 

1.  The  Gruneisen  parameter  is  defined  thermodynamically  as 

Y  =  /dP\  _  SKS 

pcp  \oT)  P  \3V/  S  pcp 

For  an  isotropic  elastic  solid,  this  becomes 


Substituting  appropriate  values  for  aluminum  gives 
=  1.97 

2.  Slater  developed  a  relation  that  was  an  extension  of  the 
Debye  Theory  for  a  three-dimensional  continuum,  and  as¬ 
sumed  Poisson's  ratio  was  independent  of  volume  and  y  is 
the  same  for  all  vibrational  frequencies: 


which  simplifies  at  zero  pressure  to 

’o  *  T  ("os  -  j)  •  2S  ■  § 

where  S  is  from  U_  -=  Cn  ♦  Su  . 

O  O  p 
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Ultrasonic  data  gives  S  =  1.40  and  y  =  2.13. 

Shuck  wave  data  gives  S  =  1.4  and  y  =  2.13. 

3.  Dugdale  and  MacDonald  give  an  expression  for  plane  wave 
velocity  vs.  volume  in  a  three-dimensional  continuum  that 
leads  to  an  expression  for  y  differing  from  Slater's  by 
a  constant  (1/3) : 


which  simplifies  at  zero  pressure  to 


Both  ultrasonic  and  shock  wave  data  give  y  =  1.80. 

4.  Anderson  and  Dienes  give  expressions  for  y  in  terms  of  the 
longitudinal  and  shear  modes  of  wave  propagation: 


Ultrasonic  data  gives  y^  =  2.31  and  y g  =  2.70. 

5.  Schuele  and  Smith  give  expressions  for  the  " low- tempera¬ 
ture"  and  "high- temperature"  Gruneisen  ratio,  where  low- 
temperature  is  the  range  in  which  it  is  assumed  the  solid 
behaves  as  a  continuum  and  high- temperature  is  the  range 
where  all  vibrational  modes  are  excited  and  y  is  constant: 
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YLT  2+£  3  (  e  yL  +  2ys) 

(£  -  W 

YHT  3  (YL  =  2yS 

Using  ultrasonic  data  and  yl  and  yg  from  above  gives 
YLT  =  2*68  and  yht  =  2*57, 

6.  Perry  calculated  the  Grueisen  parameter  from  laser  inter¬ 
ferometer  data  obtained  for  electron  beam  energy  deposi- 
* 

tion.  For  polycrystalline  aluminum,  he  obtained 
y  =  2.02  assuming  instantaneous  deposition  and  y  =  2.11 
assuming  50  nsec  deposition  time. 

(All  of  the  ultrasonics  data  used  in  this  appendix  are  for 
zero  pressure  and  20°C) . 


Perry,  F.  C.  "Thermoelastic  Response  of  Polycrystalline 
Metals  to  Relativistic  Electron  Beam  Absorption",  J.  Appl. 
Phys . ,  Vol.  41,  p.  5017-5022,  1970. 
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APPENDIX  C 

GRUNEISEN  PARAMETER,  PRESSURE  DERIVATIVE 


1.  Thermodynamic  Relation  - 


This  can  be  simplified  by  using: 


(  3P  )  T  3  (  KT  )  (  »  )  ■> 


See,  e.g.  Thurston,  R.  N.,  "Ultrasonic  Data  and  the  Thermo¬ 
dynamics  of  Solids",  IEEE  Proc. ,  Vol.  53,  p.  1320-1336,  1965. 
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to  give : 


This  can  be  modified  to: 


2.  Constant  y/V: 
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At  P  =  o,  this  reduces  to 
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A.  The  flow  stress  at  low  strain  rate  from  the  uniaxial 
stress  tests  (see  Figure  3)  is  ^  3.4  kbar,  which  gives 
=  2.3  kbar. 

3.  For  a  hugoniot  elastic  limit  )  of  5.7  kbar,  the  rela¬ 
tion  Y  =  ;  (l-2v/l-v)  gives  Y  =  2.9  kbar.  From  Equa- 

G 

tion  D-2,  =  1.9  kbar. 

C.  Since  can  also  be  taken  as  the  difference  between 
D 

~  and  P  at  a  compression  corresponding  to  the  hugoniot 
ii  H 

elastic  limit,  one  has 


~D  =  CH  '  PH  =  ;!e  ‘  Pe 


where 


=  VY2  (irf) 

te)  -  (  a.) 


Pe  =  °0CB 


For  c  =5.7  kbar:  u  =  0.0052  and  P  =3.8  kbar.  This  gives 
e  e  e 

oD  =  1.9  kbar. 


The  agreement  of  these  estimates  with  the  value  of  2.3  kbar 
from  Equation  D-l  is  considered  good,  taking  into  account 
the  assumptions  underlying  each  method  and  the  experimental 
inaccuracies  - 


Also,  it  should  be  noted  that  the  stress-axis  intercept  is 

larger  in  the  -  u  plane  .han  the  r  -  u  plane.  Assuming 
J  n  1  H  p 

Linear  clastic  and  hydrostatic  behavior  up  to  the  elastic 
limit,  then: 
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°H  "  Up  plane  " 


°D  '  °0Ue^CL  ~  CB^ 
u 


o„  -  y  plane  - 

n 


D  1+y 
y  e 


( P0C 


P0CB 


) 


Using  the  approximations  =  y  /CT  and  y  =  y  /(1+y  ), 

6  G  Li  0  0  0 

and  subtracting: 


A  a 


D 


For  a 
to  Aa 


D-l. 


5.7  kbar,  y  = 

1 . 3  kbar  from  a 


0CBUe(CL 


CB> 


0.0052  and  AaD 
H  -  Up  and  °H  ' 


=  0.S  kbar,  compared 
y  given  in  Equation 
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APPENDIX  E 

EQUATION  OF  STATE,  ULTRASONICS 


Numerous  analytical  and  empirical  relations  have  been  developed 
for  determining  equations  of  state  from  ultrasonics  data.  Seme 
of  these  are  presented  here. 


1.  The  Birch  equation  of  state  for  the  isotherm  can  be 

* 

expressed  as: 


3  T 
PT  =  2  K0 


where 


7/3  5/3 

(  /  v  2/3  ) 

(1+u)  -  (1+m) 

i-r,  (h)  -  1  } 

12-3  K 


T  1 
0T 


PT  =  1050 


7/3  5/3 

(1+  }  "  (1+u) 


2/3 


1+0.653  {(!+>)  -I 


Another  parameter  sometimes  used  with  this  equation  of  state 
is 


However,  this  term  can  usually  be  neglected,  contributing 
0.61  at  100  kbar  for  aluminum. 
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2.  The  Murnaghan  equation  of  state  is  derived  by  assuming 
the  bulk  modulus  is  linear  with  pressure  which  gives: 


Isotherm  -  P„  = 


T 

0_ 

T ' 

OT  L 


(1+u) 

r 


K 


T' 

OT 


-  1 


=  143.7 


4.87 

(l+u)  -  1 


K 

Adiabat  -  P  =  ~r 

K0S 


Ps  =  158.3 


(1  +  m) 


K 


S  ' 
OS 


-  1 


4.60 

(l+u)  -  1 


3.  Instead  of  a  linear  bulk  modulus-pressure  relation,  Keane 
assumed : 


S'.  ...  S' 

where  is  the  limiting  value  of  K  at  very  high  pressure. 


This  leads  to: 


0+  ) 
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*  S ' 

Huang  gives  the  upper  bound  for  K  as : 


KS'=  4S  -  ~  =  '  -  •=■  =  1.93 


This  gives: 


Ps  =  896 


(l+u) 


1.93 


-  1 


-  1004  In  (l+u) 


Huang,  Y.  K. ,  "Analytical  Approach  to  the  Shock  Compressibility 
of  the  18  Cubic-Lattice  Metals",  J.  Chem.  Phys.,  Vol.  53, 
p.  571-575,  1970. 
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APPENDIX  F 

WAVE  PROFILES 

(Quartz  Gage  and  Velocity  Interferometer) 


6061-T6,  6.4mm  Stock 
Velocity  ■  0. 266mm/jjsec 
Thickness  =  1. 95  mm 
Tilt  =  21ns 


] _ I _ I 

0.6  0.8  1. 

(/Li  s  ec ) 


6061-T6,  3.2  mm  Stock 
Velocity  =  0. 284mm//i  sec 
Thickness  =  3.07  mm 
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